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Abstract 

We determine effective lattice actions for the Polyakov loop using inverse 
Monte Carlo techniques. 

1 Introduction 

As there are different notions of effective actions let us start right away with the definition we will employ 
throughout this presentation: Given an action S = S[U] for some 'microscopic' degrees of freedom U 
we define an action for 'macroscopic ' degrees of freedom X = X[U] via the functional integral 

g-5eff[x] ^ j x>u6{X - X[U])e-^^^^ . (1) 

Hence we 'integrate out' U in favor of X which will guarantee that the actions S and Scs have the 
same matrix elements for the remaining degrees of freedom X. Typical examples of such effective 
actions are obtained if X corresponds to low-energy degrees of freedom, like in chiral perturbation 
theory. Alternatively, X may represent some order parameter (field) as is common in Ginzburg-Landau 
theory, for instance. 

While this is all fairly straightforward in principle one encounters difficulties in practice: in gen- 
eral, the ?7-integration cannot be done analytically. Fortunately, there is a particularly elegant way out, 
encoded in the 'effective field theory' program. There, one argues that the effective action S'eff [X] should 
have the same symmetries as the 'parent' action S[U] which suggests the ansatz 

S,s[X] = J2^kSk[X], (2) 

k 

representing a systematic expansion in symmetric operators Sk of increasing mass dimension (multiplied 
by inverse powers of the cutoff). The parameters are then fixed via 'phenomenology'. 

Alternatively, one may try to 'do the [/-integration numerically' on a lattice. However, on a lattice 
one can only calculate expectation values. So the question arises: How can one find effective actions 
from the latter? The answer is given in the next section. 



2 Inverse Monte Carlo 

Inverse Monte Carlo (IMC) avoids the functional integration in ([B by a detour consisting of three basic 
steps: first, generate configurations i = 1, . . . ,N, via standard MC procedures. Second, calculate 
the configurations {X^*) } and compute the expectation values {0[X])s = iV~^ Y^^=iO[X^''-'>]. 
Note that the X variables are distributed according to exp(— S'), i.e. the parent action. Third, in the IMC 
step proper, determine the effective couplings A^. via Schwinger-Dyson equations (SDEs) |lj||2||. The 
latter step requires some explanations. 

Generically, the target space (where the X variables 'live') will have some isometry leaving the 
metric invariant, Cx 9"'^{X) = 0, C denoting the relevant Lie derivative. This entails an invariance of 
the functional measure, VX = Y[ dX y/g{X), leading to path integral identities that can be cast into the 
SDEs 

Y,{FCxSk)sXk = {CxF)s , (3) 
k 
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where the ansatz Q has been utiUsed. The former constitute an (overdetermined) hnear system for the 
Afc which can be solved numerically. The arbitrary function if properly chosen, can be used to fine- 

tune the associated numerics. An expensive (!) check of the solution can (and should) be done by testing 
the equality of the matrix elements calculated with both the parent and daughter actions, {0)s^ff = {0)s- 

3 Polyakov-Loop Dynamics 

In f3^. W\ we have considered some particular examples of effective actions for the Polyakov loop based 
on suggestions in [5^^ JJ. 

3.1 Generalities 

At the risk of boring the experts we briefly review the physics of the Polyakov loop. This quantity is all- 
important for finite temperature Yang-Mills theory as it constitutes a gauge invariant order parameter for 
the confinement-deconfinement phase transition. The Polyakov loop is a traced holonomy (or Polyakov 
line), 

1 

[U] = —trFVa. [U] , [U] =l[Ut,a.,o. (4) 

t=l 

Introducing the Wilson coupling by = 2Nc / the theory assumes its confinement phase for j3 < j3c 
with (L) = and its deconfinement phase for (3 > f3c with nonvanishing expectation value, (L) ^ 0. The 
latter phase is characterised by a spontaneously broken centre symmetry generated by transformations 
under which L transforms nontrivially, Lx z Lx, z G Z{Nc)- 

The critical behaviour is characterised by the Svetitsky-Yaffe conjecture |8||9l which states that 
the effective theory describing finite-T SU (Nc) Yang-Mills theory in d + 1 dimensions is a Z{Nc) spin 
model in d dimensions with short-range interactions. For SU{2) this is well established on the lattice via 
comparison of critical exponents flOlfTTl . 

For gauge group SU (2) , which we consider henceforth, L is a real number between — 1 and 1 . This 
target space is somewhat nonstandard so that its isometry is not obvious. We therefore generalise to the 
group-valued variable, *p = P'^a^ € SU (2). As SU (2) = S^ its isometry is clearly an 0(4) symmetry 
generated by 'angular momenta' M^^ . Gauge invariance invariance then implies that the effective action 
can only depend on the zeroth component of *p, 5'ofr[^] = S'cff [P'^] = S'eff [L]. Restricting our identities 
io = L results in the following SDEs for the ansatz 

^^((1 - Ll)GSk,x)\k = ((1 - Ll)G,x - 3LxG) , (5) 

k 

where G^x = dG/dLx etc. Again, the functional G is useful for fine-tuning. An optimal choice is to 
use the derivatives of the effective action, G G {Si^y}, which are the operators appearing in the equation 
of motion for L. In this case the SDEs become relations between two-point functions of the Polyakov 
loop L. Again, they represent an exact, overdetermined linear system for the couplings A^. We still have 
not decided upon the concrete form of the effective action. This is somewhat of a problem, as the former 
is only mildly constrained by centre symmetry. The additional fact that L is dimensionless allows for a 
plethora of operators Sk- Certainly, a further guiding principle is needed. 

3.2 Effective Action from Character Expansion 

Recall that any (gauge invariant) function on SU{Nc) can be expanded in terms of an adapted 'Fourier' 
basis, the group characters for representation R, xr[U] = tr^ U, with U G SU{Nc)- For SU (2) we use 
the notation Xp with p = 2j and j = ^, 1, |, . . . , denoting 'color spin'. The lowest characters are the 
polynomials xi = 2L, X2 = 4L^ — 1,X3 = 8L^ — 4L, .... If we impose Z(2) symmetry and reducibility 
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Fig. 1 : The MF predictions for the VEVs of the Polyakov loop L = X\l2 (left) and the adjoint character X2 (right) as a function 
of the couplings An and A22- 

of product representations we obtain the following systematic expansion for the effective action (see also 

(xy) PI f peven 

p— qcvcn 

This action is a Z(2)-symmetric sum over higher and higher SU{2) representations including both 
nearest-neighbour (NN) hopping and potential terms. Before determining this action numerically it is 
worthwhile to study its consequences in mean-field approximation (MFA). 

3.3 Mean-Field Approximation 

If we restrict ourselves to the fundamental and adjoint representations (j) = 1, 2) and to hopping terms 
only we end up with the effective action ^ truncated such that it contains only couplings An and A22- 
This action will be denoted T2 and implies a MF potential 

p=l,2 

where z is the single-site partition function. The two MF (or gap) equations are dVuF/dxp = 0, 
p = 1,2, the solution of which yields the vacuum expectation values (VEVs) Xp of the characters. The 
latter are displayed in Fig. Q The behaviour of Z as a function of the coupling An implies the typical 
second-order phase transition (for A22 sufficiently large). Interestingly, the adjoint character X2 displays 
discontinuous behaviour as a function of A22 (for An sufficiently large). This behaviour, however, is 
located off the 'physical region' marked by the arrows. These point towards the couplings obtained by 
numerically matching to the Yang-Mills action via IMC. This is our next topic. 

3.4 Inverse Monte Carlo 

We have simulated SU{2) Yang-Mills theory on a 4 x 20^ lattice for Wilson coupling /? ranging between 
1.2 and 3.5. The critical coupling is /3c = 2.30. The associated ensembles have been matched to effective 
actions Sp with representations p < 3, hence including up to five character terms and, accordingly, five 
effective couplings. Some results are presented in Fig.|2l 

We see that both fundamental and adjoint couplings (An and A22) jump near the critical (3. The 
effect of including potential terms and higher representations (5*2, ^3) is rather mild and visible mainly 
in the broken phase {() > (3c)- The linear behaviour of An there is predicted by perturbation theory 1 13 1. 
Finally, the MF critical coupling for an Ising type effective action Si (An only) is Anc = —l/2d = 
—0.17 is off by only a few percent. 

To check the quality of our effective actions we have simulated them using the IMC values for the 
coupUngs and compared simple matrix elements namely the two-point functions (see Fig.|3ll. 
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Fig. 2: Effective couplings determined via IMC. Left: fundamental coupling An, Right: adjoint coupling \22- Different 
symbols correspond to different truncated actions. Si, ... , S3, see text. Dashed lines for small /3 are strong-coupling results. 
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Fig. 3: Two-point functions for /9 = 2.2 (left) and 13 — 3.0 (right). 



We note that in the symmetric phase (/? = 2.2) the data for different truncations lie on top of each 
other and reproduce the Yang-Mills values. In the broken phase (/3 = 3.0), however, the situation is 
different. Including higher representations leads to improvement but in particular the Yang-Mills plateau 
value (hence (L)) is not reproduced. 

To remedy this fault we have performed a brute-force calculation by also including next-to-nearest 
neighbour (NNN) interactions and a total of 14 different operators |3 |. We find that generically the effec- 
tive couplings decrease with the representation label p (as above) and the distance of the sites connected 
by the operators. A comparison of the two-point functions now yields a near-perfect match (see Fig.|4ji. 



Fig. 4: Two-point functions for /3 = 2.2 (left) and /3 — 2.4 (right) using 14 NN and NNN operators ("sim"). 

4 Summary and Outlook 

We have seen that inverse Monte Carlo based on Schwinger-Dyson equations can be a powerful method 
to numerically determine effective actions. Applying it to the Polyakov loop dynamics in SU (2) yields 
Z(2)-symmetric effective actions 5'cff[L] describing the confinement-deconfinement phase transition in 
quite some detail. The matching to Yang-Mills is reasonable if NN character interactions are used and 
becomes near-perfect if one includes NNN interactions. Obvious extensions of this work will be to 
go to SU{3) and to analyse explicit symmetry breaking terms (mimicking 'fermions'). A study of the 
continuum limit will necessitate to renormalise the lattice Polyakov loop. 
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